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Abstract 

Numerical analysis indicates that there exists an unexpected new ordered chaos for 
the bounded one- dimensional multibarrier potential. For certain values of the number 
of barriers, repeated identical forms (periods) of the wavepackets result upon passing 
through the multibarrier potential. 
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1 Introduction 

Recent studies of the chaotic aspects of different complex systems have resulted in finding 
certain conditions under which the related chaos become ordered, predicted and deterministic. 
Among these one may include the early ordered chaos found in the one-dimensional logistic 
maps [T] or the optical "chaos itinerancy" in which the onset of chaos is done in a rather 
ordered manner. Another known example of ordered chaos has been shown in [H] for the one- 
dimensional periodic potential with constant slope in which the time series are produced by 
a tent map (2]. The particle which passes through such a potential is shown ^ to be under 
the influence of an ordered chaos which is effected in an unexpected drift in the opposite 
direction to that of the average potential gradient It is also shown for the case of 
zero average gradient that the smaller becomes the widths of the potential barriers the more 
increased is the transition probability and the corresponding ordered chaos 

We note that similar results were shown |2j regarding the bounded one-dimensional multi- 
barrier potential which certainly has no average gradient and in which the width of its barriers 
is inversely proportional to the number of them. It is found analogously to [H], that the 
larger is the number of barriers along the same spatial length, which means that the smaller 
are their width, the higher becomes the transition probability of this system. Thus, since the 
bounded multibarrier system was shown in |3| to be chaotic then, according to the criterion 
in [3j for ordered chaos, we conclude that this system belongs to the last class. 

We note that although the possible presence of quantum chaos is generally expected to 
exist in systems in which the classical limit is chaotic [Sj, there are, nevertheless, quantum 
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systems which show chaotic signs without classical counterparts. This has been shown in 
for the quantum one-dimensional single barrier and in ^ for the bounded one-dimensional 
multibarrier potential. Also, the appearance of chaotic behaviour in these one-dimensional 
systems in IH E] is due to their being bounded and composed of a not large number of 
barriers. Thus, in the limit of very large number of barriers arrayed along the whole axis as, 
for example, the one-dimensional Kronig-Penney system 0, one should not expect chaotic 
effects (see also the discussion in (Hj of the localization-delocalization transition). 

Our main goal in this work is to show that changing the number or (and) the width of 
the barriers results in the appearance of new ordered chaos which is effected in the form of 
periods. These are neither periods in time nor in space but periods in the number of barriers 
A^. We note that by the phrase new ordered chaos we do not mean that the observed 
wavepackets become less chaotic and complex but that the same chaotic structure is seems 
to be repeatedly observed as the number of barriers N increases by specific values. 

We note, as mentioned, that chaoticlike effects were discussed in ^ with respect to one 
rectangular barrier and it was found that the chaotic appearance of the passing wavepacket 
as well as its correlation with the initial one critically depend upon the width of the single 
barrier. Thus, it seems appropriate to extend the discussion to the multibarrier potential 
and find the conditions under which the related chaotic effects become periodically ordered. 

We use in our discussion the well known fact that a classical (or semiclassical) wavepacket 
spreads and becomes chaotic IH [3 CHI HI] when it pass through a region along which a 
system of potential barriers (or wells) is arranged. The degree of the resulting chaos may be 
determined from the correlation |T2] between the initial and final forms of the passing particle 
(wavepacket). Thus, if this correlation turns out to be small then the initial wavepacket has 
been considerably changed and its chaotic effects have increased in the passage through the 
barriers. In the following we use the Lanczos tridiagonalization method [121 CH CH] for cal- 
culating the correlation (and, therefore, the degree of chaos) between the initial wavepacket 
and the one which emerges from the potential array. 

We show, using the obtained data of the correlation, that the passing chaotic wavepackets 
are, under cetain values of and c (which is the ratio of total interval to total width of the 
potential array), strictly periodic and predicted. That is, suppose that a specific wavepacket 
which pass through a multibarrier potential with a given A^, c and total length L assumes 
some specific form. Then it is shown that the same initial wavepacket assumes exactly 
the same form with the corresponding correlation when it pass through (A^ + nP) barriers 
arranged along the same length L and c where n is the positive numbers 1, 2, 3, ... and P are 
the periods. 

We also show that although the passing wavepackets are highly sensitive to variations of 
the ratio c, there are, nevertheless, some specific values of A^ for which the waveforms and 
the corresponding correlations are preserved for whole ranges of c. 

In Section 2 we present the bounded one-dimensional multibarrier potential and the 
formalism of the energy level statistics IHHI] used for discussing its chaotic properties. This 
method were used in |1| for demonstrating that this system is chaotic. We use here numerical 
analysis for further discussing the mentioned periodic ordered aspect of these chaotic effects. 
We show the dependence of the passing wavepackets (and their corresponding correlations 
with the given initial one) upon the number of barriers A^ and the ratio c. In Section 3 
we demonstrate the mentioned order which is effected: (1) through the remarked periods P 
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which turns out to be of two kinds; one is very frequent and is effective for a large specific 
values of and the other is rare and shows up only in two specific values of N (from the 
range 2 < < 72) and (2) through the constancy of the passing wavepackets and their 
related correlations for certain ranges of the ratio c. In Table 1 we show the correlations for 
some specific values of N and for 6 different values of the ratio c and indicate the correlations 
which are periodic. We conclude in Section 4 with a summary of the main points. 



2 The correlation between the initial and final wavepack- 
ets for the bounded one-dimensional multibarrier po- 
tential 

The bounded one-dimensional multibarrier potential discussed here is supposed to be ar- 
ranged along the x axis between the points x = — 10 and x = 10. Assuming that the number 
of barriers in the system is A^ one may introduce |4J the variables a and b which respectively 
denote the total width of the A^ barriers, where the potential V satisfies V > 0, and the total 
interval among them where V" = 0. Thus, one may realize ^ that the width of each barrier 
is jj: and the interval between any two neighbouring ones is Denoting the ratio of b 

to a by c and the total length of the system a + b hj L one may express ^ a and b in terms 
of c and L as 

a=^ b=J^ (1) 
1 + c 1 + c ^ ^ 

The possible existence of chaotic properties for any bounded one-dimensional multibarrier 
(or multiwell) potential system is usually determined by applying the energy level statistics 
IHEJ- In this method one begins from the following two-dimensional matrix equation 









S12 






^0 




S21 


S22 




B2N+1 



(2) 



where A2N+1 and B2N+1 are the amplitudes of the transmitted and refiected parts respectively 
of the passing wavapacket from the A^-th potential barrier. Aq is the transmission coefficient 
of the initial wave that approach the first barrier and is the refiected part from this barrier. 
The components Su, Su, S21, and S22 are the matrix elements of the two-dimensional S 
matrix which are related to the corresponding transfer matrix Q of the multibarrier potential 
(see, for example, Eqs (21) in |lj). The energy level statistics method [H] is used by imposing 
boundary value conditions at the remote boundaries of the system. In ^ periodic boundary 
conditions are used at the points \x\ = R, where R is much larger than the total length 
L = a + b of the system, so that one obtains 



A2N+ifiR) = Aofi-R) 
B2N+if{.-R) = B^fi^R), 



(3) 
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where f{R) and f{—R) denote the wavepackets at the points x 
Using Eqs ^ one may write Eq Q as 



A2N+1 


fiR) 
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In order to obtain a nontrivial solution for the vector 
equation 

det 



A2N+1 
Bo 



R and x = —R respectively. 



(4) 



we must solve the following 
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(5) 



and have expressed the 



f(-R) fi-R) ' 

In P5 we have used for f{R) (/(— i?)) the plane wave e*^^ (e~*'^^) 
S matrix elements in terms of the known transfer matrix elements (see Eqs (15) and (21) in 
IH). As a result of these substitutions one obtains from Eq (jHI), as in |3], a complex equation 
from which the appropriate energies which correspond to its real and imaginary parts are 
derived. Figure (8) in ^ shows the level spacing distribution of these energies in the form 
of a histogram which is clearly of the chaotic Wigner type [TT] . 

In this work we use, instead of plane waves, a semiclassical complex Gaussian wavepacket 
since this kind of wave function tends easily to be deformed and becomes chaotic upon 
passing a multibarrier (or multiwell) potential |4]. Also, the semiclassical character of the 
wavepacket enables one to simultaneously discuss, as done in the following, its momentum 
and position. Note that even in the quantum regime one may introduce the coherent state 
formalism (THl CZI which allows one to simultaneously define the expectation values 
of the conjugate variables Q and P. We note that in the numerical part of this work 
all the wavepackets (denoted 0), including the initial one, are numerically and graphically 
constructed from a given complex Gaussian packet Ppacket given by 



Ppacketi^, ^0, PO) ^o) 



. -4 

WqT!' 4 6 ""o e 



1 - 2itw^ 



(6) 



where xq is the initial value of the mean position of the packet in coordinate space and po 
and Wo are the initial momentum and width in p space. The width Wq is, actually, the initial 
uncertainty in the momentum. For an effective numerical simulation the space and time 
variables were discretized IH [151 with a resolution of dx 



Y and dt 



^ so that we obtain 

50 



dx > dt which is necessary for stabilizing and steadying the relevant numerical method 
|U HH]. For the initial xq, po and wq we choose the values of xo = —10, po = 3 and wq = \. 
In the semiclassical discussion adopted here we assume that the wavepacket is associated 
with a particle of mass m where for m we assign the value of \. Thus, as in jlj, the units we 
use for x, t and p are x = t = ^ and p = mv. That is, one may realize in this scaling 
that the velocities in ^ are related to the mentioned parameters x and t by |^ = In 
order to maintain the condition of E > V, where E is the energy of the passing wavepacket, 
we assign for the constant height of the barriers the value of = 2. 

The initial wavepacket which approach the multibarrier potential is expressed as 



0(t = 0) = Re'^ (Ppacket) + ^fn^ (Ppacket), 



(7) 
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where Re{Ppacket) and Im{Ppacket) denote the real and imaginary parts respectively of Ppacket 
from Eq ^ . The initial wavepacket of Eq ^ is shown at the left hand side Panel of Figure 
1. We note that by its definition the initial wavepacket 0(t = 0) from Eq Q spreads with 
time without having to pass through any potential (see Eq (P)). We are not interested 
here in this kind of known spreading but, especially, want to track and follow the unknown 
chaotic-like deformation of the packet due to its passage through the multibarrier potential. 
Thus, we numerically follow the time evolution of the real and imaginary parts through the 
multibarrier potential and obtain the passing wavepacket as 

(pit) = Re'^{Ppacket,v) + Im^{Ppacket,v), (8) 

where Re{Ppacket,v) and Im{Ppacket,v) denote the real and imaginary parts of Ppacket after 
passing through the multibarrier potential. 

We note that the number of the different chaotic wavepackets which evolve from the 
initial wavepacket of Eq is very large. For example, by slightly changing the ratio c or by 
adding or removing even one barrier results in a completely different wavepacket compared to 
the one which corresponds to the potential before the change. Since the passing wavepacket 
becomes chaotic there is generally no rule that controls its form or the correlation C between 
it and the initial wavepacket from Eq 0. As remarked, we show in the following section 
that for certain values of and c one may predict the forms, and therefore the related 
correlation, of the passing wavepackets. 

The right hand side Panel of Figure 1 shows how the initial wavepacket from the left 
hand side expands and become deformed at time t = 6 after passing through a ten-barrier 
potential whose ratio c is unity. The influence of changing upon the passing wavepacket 
and its corresponding correlation C is further demonstrated at the left hand side Panel of 
Figure 2 which shows the waveform obtained for the same c and t as those of the right hand 
side Panel of Figure 1 but for a 15 barrier potential. Note that by increasing A^ by 5 the 
wave packet becomes more chaotic and deformed compared to that at the right hand side of 
Figure 1. The dependence of the passing wavepacket upon c is shown at the right hand side 
Panel of Figure 2 which is drawn for the same A^ and t as those of the left hand side but for 
c = |. That is, decreasing c from unity to | causes the passing wavepacket to become much 
less complex and chaotic compared to the form at the left hand side. 

We refer in the following to Table 1 which shows the correlation between the same initial 
wavepacket of Eq ([Zj) and the passing one for 40 values of A^ and for six different values of 
c: c = 4, ^, |, 1, |, 0.25. The spatial length of the multibarrier was fixed to L = 20 and 
the time at which all the passing wavepackets were calculated is t = 6 which corresponds 
to increasing 300 times the mentioned time interval of dt = Thus, at this time the 
initial wavepacket have passed through all the barriers arranged along the fixed length of L. 
Each of the tabulated values of the correlation was numerically calculated using the Lanczos 
tridiagonalization method [121 CH CHI which yields a tridiagonal matrix the values in its 
principal diagonal are the sought-for correlations. In this method the better and accurate 
result is given by the matrix element located at the bottom of the principal diagonal. Thus, 
the larger is the tridiagonal matrix the more accurate becomes the correlation associated 
with this matrix element. This is due to the large number of numerically running the 
program which generally yields better results. We note, however, that the exact values of 
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Figure 1: The left-hand side figure shows the initial wavepacket from Eq (0) as function of 
X where x is given in units of The right hand side graph shows this wavepacket at 
time t = 6 after passing a multibarrier potential composed of 10 barriers whose ratio of total 
interval to total width is c = 1. Note how expanded and deformed the wavepacket becomes. 
Also note that the multibarrier potentials are not shown in this figure and in Figure 2 
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Figure 2: The left hand side Panel shows how the initial wavepacket from Eq changes 
at time t = 6 after passing a multibarrier potential composed of 15 barriers whose ratio is 
c = 1. Comparing this graph to the right hand side Panel of Figure 1 one may realize that 
the wavepacket becomes more deformed and chaotic by increasing N from 10 to 15 retaining 
the same ratio of c = 1. The right hand side Panel shows the wavepacket obtained at time 
t = 6 for the same 15 barrier potential as that of the left hand side Panel but with a ratio 
of c = |. The decrease in c results in a wavepacket which is less complex compared to that 
at the left hand side. 
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the correlation is not of our main concern here but especially we concentrate our attention 
upon its dependence on c and N. These dependencies may be established even from a small 
tridiagonal matrix if we use consistently the same order of it for all c and N. Thus, we use 
for all the numerical work here a third order Lanczos tridiagonalization matrix. An example 
of such a matrix is the following one which corresponds to c = |, = 10, L = 20, Xq = —10, 
Po = 3, Wo = |, t = 6 and V = 2. 



M{N = 10,c = 7/3) 



2.630821735 33.07455189 
33.07455189 1396.833502 52187.77697 
52187.77697 1669.69722 



The values of the correlation between the passing wavepacket at the time t = 6 and the 
initial one from Eq ((Zj) are tabulated along the principal diagonal. The value of 1669.69722 
at the bottom of this diagonal is, as remarked, more accurate than the two other values. The 
tridiagonal matrix is symmetric which means that the off diagonal matrix elements which are 
symmetrically located about the principal diagonal are equal. These off diagonal elements 
are the normalizing factors of the diagonal elements [T^ fT5]. 

Each one of the values tabulated in Table 1 is obtained from the bottom value of the 
principal diagonal of the corresponding tridiagonal matrix. The 40 rows in Table 1 correspond 
to the three representative ranges of = 4, . . . , 15, A^ = 31, . . . , 40 and A^ = 55, . . . , 72 and 
the six columns to the six values of the ratio c = 4, |, |, 1, |, 0.25. As seen from 
the table the correlation ranges, for the specific values given here to the related parameters 
L, wq, Xq, Po, V and t, over values which greatly differ among them. Thus, in order to be 
able to graphically plot the correlation as function of A^ (or c) one have to scale the ordinate 
axis in a log basis as done in Figure 3-4. We must remark that Figures 3-5 are constructed 
not only from the tabulated values of Table 1 but also from other values which are not 
given in this table. That is, the correlation values used for Figures 3-5 are for all A^ from 
N = 4 + n, = 1,2, . . . ,68 and for c = 4, |, |, 1, |, 0.25. From Table 1 and Figures 
3-4 one may realize that the correlation changes in a stochastic and unexpected manner even 
when adding or removing only one barrier. Also, one may see that the larger values of the 
correlation C are found at either large or small values of the ratio c and the smaller values of 
C are found at the intermediate values of c. This is shown in Figure 3 in which we compare 
at the left hand side Panel of it the correlation C as function of A^ for c = 4 (continuous 
curve), c = I (dashed curve) and c = | (dashdot curve). At the right hand side Panel of 
Figure 3 we compare the correlation C, as function of A^, for c = 4 (continuous curve), c = 1 
(dashed curve) and c = | (dashdot curve). Remembering that the ordinate is scaled in a log 
basis one may realize, for example, how large is the difference for 10 < A^ < 25 between the 
correlation for c = 4 and those obtained for c = |, c = | and c = 1. Similar differences are 
demonstrated at the left hand side Panel of Figure 4 where the correlation C, as function of 
A^, for c = 0.25 (continuous curve) is compared to those for c = | (dashed curve) and c = | 
(dashdot curve). Note again the large differences for 12 < A^ < 25 between the correlation 
C obtained for c = | and those for c = 0.25 and c = |. The similarity of the correlations for 
large and small c is demonstrated at the right hand side Panel of Figure 4 where we compare 
C for c = 4 (continuous curve) to that for c = 0.25 (dashed curve). Note, however, the large 
difference between these correlations for 28 < A^ < 36 where the C's for c = 4 are much 
larger compared to those for c = 0.25. 
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Figure 3: At the left-hand side Panel we compare the correlation C as function of N for 
c = 4 (continuous curve) to those for c = | (dashed curve) and c = | (dashdot curve). The 
right-hand side Panel compare the correlation C for c = 4 (continuous curve) to those for 
c = 1 (dashed curve) and c = | (dashdot curve). Remembering that the ordinate axis is 
scaled in a log basis one may realize that the correlations forc = 4, c=|,c=|,c=| and 
c—1 widely differ from each other. See, especially, the differences, for 10 < N < 25, among 
c — 4:, c=|, c=| and c — 1. 




Figure 4: At the left-hand side Panel we compare the correlation C as function of N for 
c = 0.25 (continuous curve) to those for c = | (dashed curve) and c = | (dashdot curve). 
Taking into account that the ordinate axis is scaled in a log basis one may realize, for 
example, the large differences for 12 < iV < 25 between the correlations C for c = | and 
those for c = 0.25 and c = |. At the right-hand side Panel we compare the correlations 
C for c = 4 (continuous curve) to those for c = 0.25 (dashed curve). From the similarity 
between the two graphs (which, however, greatly differ for 28 < < 36) one may see that 
the correlations are large for either large or small c. 
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From Table 1 and the corresponding Panels of Figures 3-4 one finds the appropriate N 
and c for constructing a bounded one-dimensional multibarrier potential from which one may 
obtain a large correlation between the initial and final wavepackets. 



3 The periods of chaos 

As realized from the former section the wavepacket which passes through the multibarrier 
potential becomes deformed and chaotic and the emerging waveforms are quite different 
even for neighbouring values of N and c. That is, the resulting waveforms (pit) and the 
corresponding correlations C depend upon c and N in such a manner that slightly changing 
either one of them results in a large change of (f){t) and C. We have, nevertheless, found 
that there exists an unexpected order among the multitude of the chaotic waveforms and 
the corresponding correlations. This order is reflected in periodicities which are clearly 
observed for speciflc values of N and c. That is, we find that exactly the same identical 
wavepackets emerge from the potential barriers when the number of the latter increases by 
specific numbers P or by any integral multiplication of them where the total spatial length 
of the system and the ratio c remain fixed. The specific periods P are found to be of two 
kinds: a large period of Pl = 140 and a small one of Ps = 28. That is, if the relevant N 
(for a specific c) is periodic then exactly the same wavepacket emerge from all the potentials 
which have {N + nP) barriers where n denote the whole numbers 1,2,... and P is either the 
large period of 140 or the smaller one of 28. Note that since 5 • 28 = 140 then any N and c 
which are characterized as being periodic with the small period of 28 are also automatically 
periodic with the larger period of 140. 

The large period of 140 is found to be very frequent and common for a large number 
of different c and N whereas the smaller period of 28 is rare. The criterion used here for 
characterizing any pair of N and c as periodic is that they have the same identical tridiagonal 
matrix for all (A^ -|- nP) where n = 0,1,2,... and P is either 140 or 28. In Table 1 we have 
denoted the multibarrier potentials which are periodic with the large period of 140 by the 
word p attached to the numerical values of the corresponding correlations C. All the other 
values in Table 1 in which the word p is absent are nonperiodic. Thus, as seen from the 
table there exists a large number of periodic multibarrier potentials which produce the same 
wavepackets and the same correlations when the number of barriers are increased by 140 
or by any integral multiplication of it. It is found (see Table 1) that the smaller is the 
number of barriers A^ the more frequent is the number of the nonperiodic potentials and 
as A^ increases the periodicity of the corresponding multibarrier potentials becomes more 
common and frequent. 

In Figure 5 we have schematically drawn for the six values of c the correlations C as 
functions of A^ from the point of view of whether they are periodic with the large period of 
140 or not. That is, each periodic C, which corresponds to some epecific A" and c, is denoted 
by a point and the absence of this point for some given A^ and c signifies that the relevant 
C is nonperiodic. Note that, as remarked, the values of A^ used in this figure are not only 
those of Table 1 but all the values of A^ = 4 + n, n = 0, 1, 2, . . . , 68. From the frequent 
occurence of the gaps in the horizontal lines of Figure 5 for small A^ and from the width of 
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these gaps one may realize that there exists a large number of nonperiodic correlations C at 
these values of N. The larger becomes the more rare and narrow these gaps become which 
means that the number of the periodic correlations C increases for all values of c. For very 
high values of N (not shown in Figure 5) the continuous linear sections become very long 
for all c which means that all the wavepackets as well as their corresponding correlations are 
periodic with the large period of 140. 

Regarding the smaller period of 28 we have found that it exists for the two pairs of 
(c = i, N = 29) and (c = 0.25, N = 28). That is, the same identical tridiagonal matrix, 
which implies the same emerging wavepacket and correlation, is obtained for the c — A case 
for all potentials (arranged along the same fixed length of L = 20) which have (29 + n ■ 28) 
barriers where n ~ 0, 1, 2, .... Likewise, for the c — 0.25 case one obtains the same 
tridiagonal matrix, wavepacket and correlation (which are not the same as those of the 
formerly discussed c = 4 case) for all potentials which have (28 + n ■ 28) barriers arrayed 
along the same fixed length of L = 20 where n — 0, 1, 2, .... 

Another kind of ordered regularity which we have found among the multitude of all the 
chaotic wavepackets is related to some specific tridiagonal matrices (and, therefore, to their 
corresponding wavepackets and correlations) which remain constant even when the ratio c 
changes. That is, although the general behaviour is the unexpected change and deformation 
of the passing wavepacket when c changes even slightly there exist, nevertheless, specific 
values of N which are characterized as related to correlations and wavepackets which retain 
their forms even when c is changed. For example, for N — 70 we have found that the same 
tridiagonal matrix (which means the same wavepacket and correlation) remains constant 
for 1 > c > 0.25. The same situation is also encountered for N = 71 where this time the 
constancy of the matrix, wavepacket and correlation are retained for the larger range of 
4 > c > 1. The relevant tridiagonal matrix for the last case is 



M{N = 71,1 < c < 4) 



0.06426037493 6.899198452 
6.899198452 68.29688359 562.8969638 
562.8969638 20.28518401 



The same situation is again encountered for other values of N for which one finds constant 
different matrices . These iV's and the corresponding ranges of c over which the emerging 
wavepackets (and the appropriate correlations C) retain their forms are; at N — A7 for 
I > c > |, at = 24 for 4 > c > I and at = 18, AT = 24, AT = 53, and A^ = 123 for 
1 > c > |. Note that all these values of A^ are also characterized as being periodic with the 
large period of 140. Thus, one may realize that the constancy of the relevant tridiagonal 
matrices (and the corresponding wavepackets and correlations) are retained not only for 
these specific A'^ but also for all the other A^'s obtained by increasing them by 140 or by 
any integral multiplication of it. In other words, this behaviour of constant waveforms (and 
correlations) when c changes is strongly related to the previously mentioned behaviour of 
periodic waveforms (and correlations) when the number of barriers, for these specific A'^ and 
c, increases by 140 or by any integral multiplication of it. 
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Figure 5: From this figure one may compare and find if the 
correlations C, for the six values of c = 0.25, |, 1, |, |, 4, 
have the large period of P = 140. The six values of c label the 
6 horizontal lines. A periodic C, with given c and A^, is denoted 
by a point and the absence of this point signifies a nonperiodic 
C. Thus, as seen from the figure the larger is N in the abcissa 
the longer become the linear continuous sections for all c which 
denote that the number of periodic C becomes large. The wide 
gaps at these lines for small signify that there exists a large 
number of nonperiodic C at these values of A". 
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4 Concluding Remarks 

We have discussed the chaotic deformed wavepackets which come out of a bounded one- 
dimensional multibarrier potential and study the correlation of these wavepackets with the 
initial one. It has been shown, using the Lanczos tridiagonalization method, that there 
exists an unexpected order and regularity among the multitude of all the possible chaotic 
wavepackets which come out of this system. This order is characterized by the existence of 
two periods through which one may obtain the same wavepackets and correlations when the 
number of barriers increase, for some specific A^" and c, by either 140 or 28 or by any integral 
multiplication of them. The more common and frequent period is that of 140 whereas the 
smaller one of 28 is rare. Any wavepacket and its corresponding correlation which is periodic 
with the small period is also automatically periodic with the larger one. 

The correlation C, as function of either N or (and) c, between the passing wavepackets 
and the initial one is stochastic and discontinuous as may be realized from Table 1 and 
Figures 3-4. One may see from these figures and from Table 1 that the larger values of C 
are obtained for c ~ 4 or c 0.25. 

Another ordered behaviour that we have found is related to the constancy of the wavepack- 
ets and the corresponding correlations for specific A^ and for whole ranges of c that may be 
as large as Ac = 3. All these A'"'s are also characterized as being periodic with the large 
period of 140. 

In summary, one may see that the chaos demonstrated by the bounded one-dimensional 
multibarrier potential is an ordered and periodic phenomenon especially for large A'^. 
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Table 1: The table shows the correlations C between the passing wavepacket at time t = 6 
and the initial one from Eq 0. The rows correspond to the three ranges of = 4 + n, n = 
0, 1, . . . 11, = 31 + ra, n = 0,1, . . .9, and N = 55 + n, n = 0,1,... 17. The columns 
correspond to the ratio c for c = 4, |, |, 1, |, 0.25. Any which have correlation C with 
the word p attached to it is periodic with the large period of 140. That is, any such N have 
exactly the same value of C for all multibarrier potentials which have (A^ + n ■ 140) barriers 
where n = 1, 2, 3 . . .. 
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